An exact property of the nonequilibrium photon Green function for bounded media 
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The nonequilibrium photon Green function for a bounded medium surrounded by vacuum is 
analyzed on the basis of the Dyson equation. As its components, the field-field fluctuations as well 
as the spectral function split up into parts related to medium and vacuum. Particularly, it is shown 
that the vacuum-induced fluctuations describe propagation of arbitrary, even nonclassical light in 
terms of solutions of the classical wave propagation problem. The results generalize previously 
obtained ones for steadily excited media in slab geometry. 

PACS numbers: 



I. INTRODUCTION 

Recently, we were able to present a general, exact law 
for the electromagnetic energy flow between a medium 
and the surrounding vacuum The dynamical and spa- 
tially inhomogeneous coupled system of matter and light 
was described in a quantum-kinetically consistent way 
by photon Green functions, i.e., Green functions (GFs) 
defined in terms of the electromagnetic vector potential, 
and the nonequilibrium Keldysh technique (see Ref . for 
a short introduction to these topics). 

Many exact relations could be derived, which can help, 
e.g., to describe emission and lasing in excited semicon- 
ductors, spectral signatures of quantum condensation, or 
propagation of nonclassical radiation.- However, the va- 
lidity of the obtained relations could only be shown as- 
suming the steady state and a highly symmetric geome- 
try, namely a medium slab with finite thickness but in- 
finitely extended in the other dimensions. 

Before, this technique was primarily applied to in- 
finitely extended (bulk) systems, and thermodynamic re- 
lations as well as kinetical equations have been estab- 
lished for this case. However, little attention was devoted 
to the fact that real-world matter is never infinitely ex- 
tended but instead exchanges energy with the surround- 
ing environment. This aspect applies the more if optical 
emission and absorption of bounded media are consid- 
ered. 

All the results mentioned above are closely related to 
a surprizing property of the photon GF, namely its strict 
splitting into vacuum- and medium-related contributions 
which enter the physics of emission and absorption in a 
completely different way. Whereas the medium-related 
ones are well-known from the optical theorem for bulk 
media, the vacuum-related contributions result from a 
subtle cancellation between infinitesimally vanishing and 
diverging factors. This has been shown in Refs.QJ^ under 
the mentioned restrictions, but from physical intuition 
it is clear that this behavior cannot depend entirely on 
some special assumptions but should instead reflect an 
universal property. 

This applies as well to another, rather counter-intuitive 
property of the photon GF shown in the mentioned work, 
namely that certain quantum- mechanical effects (e.g., 



emission) are completely determined by optical proper- 
ties of the media which originally describe propagation 
of classical light (e.g., absorption). 

In this paper, I will analyze the photon GF further 
and show that these results can be generalized to some 
extent to complete universality. Indeed, the splitting will 
be shown to be a universal and exact property of the 
photon GF for bounded media, which, to the best of my 
knowledge, has not been reported before. Furthermore, it 
follows that the propagation of field fluctuations through 
arbitrarily dispersive and absorbing media can be com- 
pletely described in terms of the solution of the classical 
wave propagation problem. 

The exact proof will be given in the following. In con- 
trast to previous approaches, it starts with the integral 
form of the Dyson equation for the Keldysh photon GF, 
because only in this representation the contributions in- 
duced by the vacuum appear quite naturally. With the 
help of the Langreth theorem^ expressions for the field- 
field fluctuations can be extracted which already exhibit 
the desired formal structure (splitting). 

Further investigation of the vacuum-related contribu- 
tions requires a re-inspection of the vacuum case. Note 
that the term "vacuum" in this paper will denote the ab- 
sence of a medium, where the polarization tends to zero. 
However, electromagnetic radiation may be present in 
this vacuum or free space, and it is thus not necessar- 
ily an electromagnetic vacuum. Consequently, the term 
"vacuum fluctuations" will denote quite generally all pos- 
sible fluctuations of the electromagnetic field in the free 
space. Those appearing already in the ground state, i.e., 
without any external preparation of the free space, will be 
called "spontaneous" fluctuations, while fluctuations in- 
duced, e.g., by preparation of the free space as a (photon 
number) Fock-state or even thermal (hohlraum) radia- 
tion, are called "stimulated" ones. 

First it is demonstrated that already the spectral func- 
tion of the vacuum can be written as a product of an 
infinitesimally vanishing term with a just equally strong 
diverging integral. Using this in the formal structure 
for the field-field fluctuations proves that the vacuum- 
induced contribution to the full spectral function can be 
given in terms of the retarded (advanced) GF governing 
the propagation of classical light. 
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Then the vacuum-induced field fluctuations are stud- 
ied starting with a normal mode/plane wave expansion 
for the free photon GF. With the aid of it, propagation 
of arbitrarily prepared light can be traced back to the 
solution of the classical wave propagation problem. 



II. DYSON EQUATION 

The following summary of the most important GF def- 
initions follows closely the more extensive presentation in 
Ref. i The photon GF in a many-particle system cou- 
pled to electromagnetic fields is defined in terms of the 
vector potential operator A. The effective vector po- 
tential in Coulomb gauge on the double time contour, 
A(l_) = (A(l_)) c , where 1 is the parameter tuple (ri,^) 
and t denotes a time on the contour, obeys Maxwell's 
potential equation in its usual form, 

-Mo [jind(r,£) +Jext(r,£)] , (1) 

where the transverse vector fields j; n( j = (j) c and j ext de- 
note the medium-induced and externally induced current 
densities, respectively. Note that the latter is a c num- 
ber function. Additionally, Ao and A ex t are introduced 
as the solutions of the free homogeneous and inhomoge- 
neous wave equations: 



A -^l Anlr. /■)=!.. 



i d 2 



(2) 



A ~ # gj5 ) A ^t (r, t) = -^ojoxt (r , t) (3) 



The photon GF D(l, 2) is defined as the functional 
derivative 



MO <5jext,fc(2) 

~{(^^) c -{^)(^&i}}, (4) 



Afc(l»2)=- 

1 i 

Mo 



which contains the above-mentioned operator correla- 
tions {i,k denote vector components). Accordingly, the 
polarization function P ("photon self-energy") is given 
as 



P ife (l, 2) = -Mo 



Sjindd (1) 

6A k (2) ■ 



(5) 



From these two definitions, one can obtain the Dyson 
equation, which reads in its integral form 



Afc(I,2) 



Afc,o(I,2) + £>ij,oa,3)Pji(2,4)Afc(4,2), (6) 



where Dq is the free, or vacuum, photon GF, whose in- 
verse is 

^, 1 o(L2)=(A ri -lg) ( 5(l-2). (7) 

Here and in what follows, the sum convention is ap- 
plied. GFs on the double time contour contain four 
different physical functions ( "Keldysh components" ) ac- 
cording to the possible operator time orders on the con- 
tour. Only two of them are independent, and one usually 
chooses the "greater /less" components and defines 
the retarded and advanced GFs 

Df (1,2) = 9(tr -t 2 ) {D>(1,2) - D<(1,2)} , (8) 
{1,2) =Df {2,1), (9) 

which govern the classical wave propagation problem 2 
according to (in short-hand notation) 



D *et,-i A = l D ret,-l _ p ret . A = 



(10) 



This can be rewritten as the Lippmann-Schwinger equa- 
tion 



or finally as 



A = A + A ext + £>fP ret A 



ret,-l 



(11) 



(12) 



where e™* is the transverse dielectric tensor, 



SA C 



Jjext SA C 



vet _ SA cyLt _ t t 



(14) 



According to the Langreth rules^ any product of 
Keldysh GFs F lk {l,3)G kj (3,2) has the "greater/less" 
Keldysh components 

[FG)%{\,2) = 

F%\1, 3)G|(3, 2) + F*(l, 3)G$ V (3, 2) (15) 

and the following product rule applies for the retarded 
and advanced GFs, : (FG) rct / adv = fct/advgrot/adv^ 

With the help of these rules, one obtains after some 
rearrangements as formal solutions of the Dyson equation 
([6|) the field- field fluctuations (in short-hand notation): 



^ ~ v i 

£>! = D' ct P^D aAv , 



(16) 
(17) 
(18) 



Eq. (|16|) generalizes the optical theorem (also referred 
to as dissipation-fluctuation theorem in particle kinetics) 
to bounded media. The field-field fluctuations D< split 
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up into a well-known contribution proportional to the 
medium polarization function, the medium-induced con- 
tribution D^, and, additionally, the vacuum-induced con- 
tribution Dv , which can be related to the vacuum or free 
space fluctuations Dq. 

These two kinds of field-field fluctuations were shown 
to enter the physics of absorption and emission in a com- 
pletely different wayJ^ In the steady-state slab geome- 
try, namely, the vacuum-induced contributions alone de- 
termine absorption and emission, while medium-induced 
contributions completely cancel out. Here, no restrictions 
or assumptions were made, so this splitting appears nat- 
urally as an universal property of the photon GF, and 
is valid for arbitrarily nonstationary excited media with 
any geometry. However, the general physical role of the 
contributions may not trivially be deduced from the spe- 
cial case of steady-state slab geometry. 

III. SPECTRAL FUNCTION 

The spectral function is introduced as 

A, = Df Df = D>~D<. (19) 

Here the first identity should be regarded as the defini- 
tion of D, which is unique due to the radiation condition 
for £) rct / adv (outgoing/incoming radiation). The second 
identity follows as an identity between the various com- 
ponents of the Keldysh GF and is crucial in the following. 

The spectral function equally splits into medium- and 
vacuum-related contributions D m and D v according to 
Eq. (|16[) . In the latter contribution, the spectral function 
of the vacuum enters, 

D = - Df v = D> - D< , (20) 

which is given exclusively by vacuum functions. These 
are spatially and temporally homogeneous functions and 
thus depend on their difference variables only, e.g., 
£>Q Ct ^ adv (l, 2) = Z?Q Ct / adv (l — 2), so that one may Fourier 
transform (ri — r2, t\ — t-i) to (q, uj). Now, all the integral 
relations reduce to algebraic ones. 

Also, due to isotropy in these functions, their tensorial 
character reduces to the product of a scalar function with 
the transverse unity tensor, e.g., 

D^.(q,L,)=t !3 (qX(q, W ), (21) 

= ^2 eAq, 4 e Aq j = - (22) 

A=l,2 9 

where e Aq denote the transverse polarization vectors. 
The function Dg et in Fourier space is given by 

D^iq, u>) = — — -| 5-j = D adv (q, w)* , (23) 

[ui + icy — c z q z 

where e — > ensures the causal structure. With this 
relation inserted, the spectral function of photons in the 



vacuum case, Eq. (|20[) , can be written in the well-known 
form 

7TC 2 

A)(q, to) = -r— [S(uj - cq) + 5{uj + cq)} 

nc / ( 24 ) 

= — \5(u) - cq) - 8(to + cq)] , 

which follows, e.g., by applying Dirac's identity 
(Sokhotsky-Weierstrass theorem) to Eq. (f2"3"| . Moreover 
and non-trivially, this can shown to be identical to 

D (q, to) = - ie ^f(q, c)D adv (q, co), (25) 

a product of the infinitesimally vanishing factor e — > 
and the correspondingly diverging term 

^ (q ,^( q ,^ ( ^ c2 ;; 2+4 ^ 

4 

= — r [6(u-cq) + 5(u + cq)]. (26) 

Now in the original (space-time) domain Eq. (|25|) can 
be written 

D (l, 2) = ^D^(l, 3) A D§ dv (3, 2) , (27) 

which together with (fT3| used in Eq. (fT8|) yields 

D v (l, 2) = ^D rct (l, 3) J- D adv (3, 2) . (28) 
<r ot 3 

Just as Do, D v shows up as a product of e — > times an 
improper integral diverging with just the same strength 
oc 1/e. However, this integral now consists of the full 
retarded/advanced GFs instead of the vacuum ones. 

IV. VACUUM-INDUCED FLUCTUATIONS 

Whereas the spectral function of the vacuum is 
fixed, its kinetical state (or population of modes) is still 
to be specified corresponding to the (external) prepara- 
tion of the outside region. One may thus assume, quite 
generally, an arbitrary field-field fluctuation D^ stim (l, 2) 
in the vacuum or free space. It is, hence, externally 
stimulated by preparation and will add to a spontaneous 
contribution from the ground state fluctuations of the 
electromagnetic vacuum, D^^, such that the complete 
fluctuations in the vacuum can be written as 

£>£(1, 2) = £>J sp (l - 2) + £>| 8tim (l,2). (29) 

It is insightful to apply in D^-(l, 2) a normal mode ex- 
pansion for the vector potential operator 5 following the 
definition (|2|), e.g., for a vacuum Fock-state with the pho- 
ton population n\q = = — 1. One then finds the 
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tensor 
where 



(1) ®F> q (2)* 



n? q F^(l)®F< q (2r], (30) 



Ffqi 1 ) = e Aq ex P [i(qi"i =F cqh)] 



(31) 



describe classical plane waves with polarization e> q and 
wave vector q. 

From the Fourier domain representation 



££(q,w) 



7TC V 



[e\d ® e Aq ] 



x - «?) + ™Aq<H w + ( 32 ) 

one easily obtains the ground state contribution (n = 0), 



1? 



: 5(w =F eg). 



(33) 



Note also, that according to Eq. (f2TJ)) . the vacuum spec- 
tral function (|24|) is, as it should be, always conserved in- 
dependently of the choice of D> stim (l, 2) = D< stim (l, 2). 

Using Eq. (|18|) , the vacuum fluctuations appear renor- 
malized due to the presence of a bounded medium ac- 
cording to 



0(1,2) 



Y — 



^ q A> q (l)®A> q (2)* 



'Aq ^Aq 

where the effective fields 

5= ret,-l 



nlA<„(l) ®A< q (2)*], (34) 



A^ q (l) = 4 C ^ i (l I 2)F< q (2) 



(35) 



describe propagation (i.e., reflection, absorption, and 
transmission) of a classical plane wave F^- q in the pres- 
ence of a bounded medium. 

A special case is the preparation of the vacuum to ther- 
mal equilibrium or hohlraum radiation, for which due to 
the Kubo-Martin-Schwinger condition n\q becomes the 
Bose distribution^ 



"Aq 



exp 



Hcq 



- 1 



(36) 



In Ref.0, Eq. 37, a nonequilibrium distribution for the 
vacuum- induced photons, n(qii,a;), was introduced. It is 
closely related to the photon population n introduced 
above. In fact, in the case of pure vacuum (P(l, 2) = 0) 
and for cylindrical coordinates q = ((fa,q||), any depen- 
dence on q may be rewritten as one on (qy , uj) due to the 
fixing of uj = ±cq through the <5-functions in Eq. 



since q 2 — lu 2 /c 2 
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V. CONCLUSION 

In the analysis of the photon Green function for 
bounded media, the splitting of the field-field fluctua- 
tions into contributions induced by the medium or the 
vacuum appears as a universal and exact property, as it 
was to be expected and is now demonstrated for the first 
time. Despite of their different physical role, the defi- 
nition of both contributions, Eqs. (fl"7|) and (fl"8]) . shows 
a common pattern in that the light propagation due to 
an arbitrary kinetic state (of the medium, P^, or of the 
vacuum, Dq) is fully determined by retarded/advanced 
GFs. These solve the classical wave propagation prob- 
lem, i.e., reflection, transmission and absorption, and are 
given by the optical properties of the medium. Thus, 
propagation of arbitrarily prepared (even non-classical) 
light can be traced back to the classical wave propagation 
problem. 

The medium-induced fluctuations D m are determined 
by D lct . Although an analytic expression for D rct could 
be given in Refs. [l|@ with restriction to steady-state slab 
geometry and some limitation to the spatial parameters 
(sources outside the medium), this seems to be out of 
reach for the general case, i.e., arbitrary geometry in non- 
equilibrium. This makes D m hardly accessible to any 
further analysis. In contrast, the transversal dielectric 
function, which determines the vacuum-induced fluctua- 
tions, is given in Eq. (O by D 1 ^ (Eq. (J23j>) and P rct . 
It is, hence, far more accessible, since P lct trivially re- 
lates to the optical susceptibility, for which analytic or 
numerical approximations can be used. 

Thus, the paper concentrates further on the vaccuum- 
induced contribution. The universality of an interesting 
identity for the vacuum spectral function, Eq. (|25[1 . was 
ensured. It translates directly to the spectral function 
of the vacuum-induced contribution, Eq. (|28[) . which, in 
this form, plays a crucial role in the derivation of the 
nonequilibrium energy flow law presented in Refs. [l]|2. 

Last, it is shown how ground state fluctuations D^ sp 
and externally stimulated light enter in this framework. 
A normal mode/plane wave expansion for the vector po- 
tential yields exact relations for D^ sp , and D^- in 
terms of classical plane waves and photon population 
numbers. 

It would be interesting to compare these relations to an 
expansion into non-classical states, e.g., squeezed light, 
and analyze the propagation of this light on this footing. 
Especially, it should be possible to replace the input- 
output formalism, a workhorse in quantum optics,— >^ by 
exact relations that account for the spatial inhomogene- 
ity inherent to bounded media problems. 
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